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■ ABSTRACT A calculation of the chiral anomaly on a finite lattice without 

OO i fermion doubling is presented . The lattice gauge field is defined in the spirit 
of noncommutative geometry. Standard formulas for the continuum anomaly are 
^ ■ obtained as a limit. 
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a^ ■ 1. THE FINITE LATTICE EFFECTIVE ACTION 



Chiral fermions have been the Achilles' heel of lattice approximations in QCD. In 
standard lattice method one uses a nearest neighbor approximation to the derivative 



^1 of a fermion field, the gauge field is represented by the Wilson link variables for 
' nearest neighbor lattice points. In this approach one cannot avoid an unphysical 
doubling of fermion degrees of freedom. In the present paper I want to propose 
an alternative lattice approximation which uses not only the nearest neighbors but 
all points on the finite lattice in the construction of the (chiral) Dirac operator. 
This may sound a bit clumsy, but actually we shall see that it is quite natural, and 
what is most inportant, it leads to the correct continuum limit without a fermion 
doubling. 

We study chiral fermions in euclidean space in 2n dimensions. In the continuum 
the fermions are functions in M^"^ taking values in ® C^. The euclidean 7- 
matrices operate on the first factor and the compact gauge group G on the second 
factor. The 7 matrices satisfy 
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The matrix 72n+i = (— ^)"^7i72 • • • 72n anticommutes with the other 7 matrices. 
We define the chiral projections P± = |(1 ± 72n+i) and consider the chiral Dirac 
operator 

(2) Da = 7m(-^^m + A^P+) =i> + 4P+ 

coupled to a vector potential JjL = ^^A^ with functions on 'M?'^ taking values in 
the Lie algebra g of the gauge group G. 

On a finite lattice L C M^"' consisting of points x = m • A, where the lattice 
spacing A is some fixed positive real number and m — (toi, TO2, • • • , TO2n) is a 
vector with integer components = —q, —q + 1, . . . ,q, any function is a linear 
combination of the Fourier modes fp{x) = (2qV" ^^^^ with = -^^k and k = 
~Qi ~Q + ^1 ■ ■ ■ iQ- Thus the vector space of massless fermions on the lattice L has 
dimension K = 2^ ■ N ■ {2q)^'^ . Note that at the end points, corresponding to A; = ±q, 
in the momentum space lattice the Fourier modes are identical; we have periodicity 
with period 2q in the momentum index k. We shall later break this periodicity. 

When studying the finite volume continuum limit one lets q ^ 00 and A — > 
such that £ = 2qA is fixed. However, we shall not insist on this but we keep the 
option open for an infinite volume limit. 

Now our approach departs from the usual setting. Normally, the Dirac operator 
is written using the principle that partial derivatives on a lattice are certain finite 
difference operators involving nearest neighbor points. Here we define the free 
Dirac operator through the Fourier transform Dq i— > ^iu,Pfj,- The vector potential 
A = A^'^^, thought of as an operator in the one-particle Hilbert space, is defined 
as a. K X K matrix with matrix elements 

(3) {lixA^)sp,a'p' = - P')lix)88' 

where ss' denote both the spin and gauge indices and the ambiguity in the Fourier 
transform A (when restricting the continuum potential to the finite lattice) is re- 
solved by requiring that both A;, k' = —q, —q + l,...,qioTPn = k ■ Note 
that in general the Fourier transform is not periodic and therefore the matrix ele- 
ments Agp^s'p' for p,p' = i^^Q are all independent. The number of lattice sites is 

/o„ I 1 ^2n ; J- ] ^'o„^27^ ] u ; ;™„1,, 1 ] 
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lattice operator A (in momentum space) as a, K xK submatrix of the corresponding 
continuum operator by restricting the momentum indices to the given range. 

We shall generalize the above setting in the spirit of noncommutative geometry, 
[C] . The reason for doing it here is purely technical: We do not need to think all the 
time whether the matrix representing a vector potential really represents a multipli- 
cation by a Lie algebra valued function. A vector potential in the generalized sense 
is any hermitean K x K matrix of the form ^^7^ where the A^'s are hermitean 
matrices acting in c^xC^Q+'^ f" ^ A gauge transformation is then an unitary K x K 
matrix g which commutes with the ^fjjs. The action of on a potential A is given 
as 

(4) 4' = 9-'49 + 9-'\^.9]- 
The curvature form is 

(5) F^,^ = [p^, A^] - [p^, A^] + [A^, A^] . 

It transforms as usual, F^^ i— > g~^F^^j^g. 
We define the lattice effective action 

(6) S,ff{A) = -logZ(^), Z{A) = detreni:^ + 4P+) 

where detj-en is a renormalized determinant. Even on a finite lattice we want to 
use certain renormalized determinants in order that in the limit q,l/A — > oo the 
effective action remains finite and leads to the continuum effective action. 

In continuum the first renormalization ('vacuum subtraction') is to replace the 
determinant det(/i + 4-P+) by det(^(^ + 4+P+)) = det(l + T), where T is an 
operator of order —1 in momenta. The term ie, for a real nonzero parameter e, is 
introduced as an infrared regularization. From now on any in the denominator 
stands for ^ + ie. Using the formula log det = tr log and the expansion of the 
logarithm 

(7) tr log(l + T) = tr(T -^T^ + ^T^..) 

one localizes the potentially diverging terms as the traces of for /c = 1, 2, . . . , 2n. 
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that the potential A vanishes more rapidly than Ixj"^" at infinity, the trace of 
these terms is finite. In terms of a momentum space cut-off A one can write an 
asymptotic expansion 

(8) trA(T - + . . . - i-T^-) = Yl ^'^i + log(A/Ao)a/o5- 

i<io 

The scale fixing constant Aq is determined by the physical requirement that the 
strength and location of the pole at g = of the boson propagator is not affected 
by the loop diagrams. The renormalized trace is defined as the coefficient qq in the 
asymptotic expansion. 

It follows that we can define the effective action as proposed by Seller [Se] , 

— logdetA;(l + T) + a finite number of Feynman diagrams 

where 

(9) logdetfe(l + T) = tr + Wt^'^' + ' ' ' ) 

with k = 2n + 1. The finite set of (renormalized) Feynman diagrams comes from 
the renormalized traces of T'' for k < 2n. For example, when 2n = 4 the only 
diverging diagrams are the vacuum polarization terms with at most four external 
gauge boson lines. 

In the lattice the modified Fredholm determinant detfc(l + T) is defined exactly 
as in the continuum case. The renormalized determinant is 

(10) det,en(l + yP+) = det2n+l(l + Up+) ' 6™^^+- + ^^^") 

where the renormalized trace TR is defined as follows. 

Case of 2n = 2. The continuum limit of trA(^^+) vanishes by a simple par- 
ity argument. Thus no renormalization is needed for this term. The next term 
trA(T^) is potentially logarithmically diverging. However, by the trace properties 
of products of 7 matrices and parity this term is actually of order —3 and gives 
a finite trace. Thus the effective action is completely determined by det3(l + T) 
and the finite 1-loop diagram, TR(T — \T'^) = — ^tr^-T^, where tr^ stands for the 
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indices, then integrates over momentum variables with \p\ < A, followed by the 
limit A — oo and the integration over x. 

Case of 2n = 4. The first term tr T vanishes as in the previous case. The next is 
the 1-loop diagram trT^ with two external boson lines. In the continuum we must 
compute traces of operators which are composed of products of Green's functions 
(the operators l/{'jl> + ie)) and and of smooth functions 4- They are examples of 
pseudodifferential operators. The algebraic manipulations involving PSDO's are 
most conveniently performed using the symbol calculus. First let us recall the 
basic rule of symbol calculus: A pseudodifferential operator is represented by a 
smooth function of coordinates and momenta, its symbol. The symbol a * 6 for a 
product of operators is computed from the symbols o, 6 of the factors as 

V J\ ^fJ ^ ni\n2\ . . .nd\ {dpiY^ . . .{dpdY" {dxiY^ . . .{dxdY"' 

Now the first term in the expansion of the PSDO (^^P+)^ leads to both quadratic 
and logarithmic divergencies which are given by the integral 

^ U^x I d^p ^ ^ tiA^ 



{2tiY j J ^(p2 + e2)2 

The lattice version of this is simply 



In addition, there is are logarithmic divergencies arising from the next terms in the 
expansion of as well as contributions from the higher order terms which involve 
the traces tr(T'^) for /c = 3,4. In the continuum case, by a standard Feynman 
integral calculation carried out in [SAB J] , one obtains {(3 is a numerical constant) 
as the total logarithmic divergence 



/51og(A/Ao) j tr(F2 + F*F) 



where Aq is a renormalization constant and *F is the dual of the field tensor. 
Actually, in our case the second term involving F * F vanishes because we are 
not considering instanton backgrounds, the vector potential is globally defined and 
vanishes at |x| — > oo. A lattice version of this is 
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modulo finite terms in the continuum fimit. Thus for 2n = 4 the effective ac- 
tion is given by Z{A) = det5(l + T) x e^^"^\ where (3{A) is the ordinary trace 
X]fc=i fc ^^{^4-+) iT^'inus the diverging terms discussed above. 

General case. In 2n dimensions there is a finite number of both polynomially 
and logarithmicaUy diverging terms. A derivative in momentum space makes the 
diagram better converging. A derivative in momentum space is associated with a 
differentiation in x space of one of the A's in the expansion of the trace. Since 
for diverging diagrams we can have only a finite number of differentiations, the 
coefficient of A'^ or of logA will be a finite differential polynomial in A. The lattice 
renormalization is obtained by subtracting these diverging terms from the naive 
effective action in such a way that the partial derivatives are replaced by the 
multiplication operators ip/j,. 

2. THE ANOMALY 

Next we shall compute the gauge variation of the (lattice) effective action. For 
that we need some properties of the generalized Fredholm determinants. According 
to our definition, [S], 

(11) detfe(l + X) = det[(l + X)e^'=W], 
where 

(12) Pk{X) = J2^—^X\ 

i=l 

If all the traces of powers of X are finite we can write 

(13) detfe(l + X)= det(l + X) ■ e^'Mx)^ 
The generalized determinants have a multiplicative anomaly, 

(14) detfc[(l + X){1 + Y)] = detfc(l + X) ■ detfc(l + Y) ■ e^^^^^^\ 
where 
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The first nonzero multiplicative anomaly is 72(X, y) = —tr{XY), the next is 
73(X, Y) = tr{X^Y + Y'^X + l{XYf). We shall also need the derivative Hf, = 7^^^ 
with respect to the first argument at {X,Y) = (0,F). Its value for a variation 5X 
is easily computed to be 

(16) Hk{Y){5X) = tr{Y''-^5X). 

The gauge variation of Seff{A) can be computed as follows. Modulo a variation 
of a finite polynomial in A, the anomaly is given by the gauge variation of the 
functional Zk{A) = detfe(l + ^4+)^ 4+ = 4^+- Denote g- = gP+ + P- and 
^+ = gPj^ + P- . We obtain 

Zk{Aa) = Zk{g-^Ag + g-^dg) 

= detfc(^(7-(|^ + 4+)^/;') 

= deik{gl^^g-i){l + ^4+)) 

(17) = deikigA-i)) ■ detfe(l + \4+) ■ ^^^i9-''3[MP+4^+)^ 

V V 

For infinitesimal gauge variations we have 

(18) ^{X-A) = bxSk{A) = ^5,(A(t))|,=o = (i?i7fe(0,^4+))(^[/^,^]^+), 

where D\ is the derivative with respect to the first argument and Ait) = A^^^^ for 
^(t) = e*^.By (16), 

(19) a;(X;A) = tr(lb^,X](l4+)'=-i). 

Example 2n = 2 Because of the accidental property of the two dimensional 
effective action, (^4^+)'^ is conditionally convergent, we may use the determinant 
det2 (instead of the more complicated dets). The anomaly is then computed from 
the multiplicative anomaly 72. Both on the lattice and the continuum the anomaly 
is 

(20) io{A-X) = tr (-l-y,,X]P+-^4P+] . 



]i> + ie ' ]l> + ie 
In the continuum a simple computation leads to 

(21) uj(A-X) = — /" (fxtTA„,d^Xe'''' + ^ ! (fxiiA,,d,,X. 
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The second term is a trivial cocycle, it is the gauge variation of the local functional 

47r J 

The first term is the nontrivial part of the anomaly. On the lattice the local 
formula does not make sense, but we still have the trace formula (20), which in 
the continuum limit leads to an integration in momentum and configuration space, 
giving the formula (21). 

The general case. The nontrivial part of the anomaly comes from the difference 

(22) OkiA) = log det(l + T) - logdetfe(l + T). 

This follows form the fact that the unrenormalized determinant is perfectly gauge 
symmetry; the gauge symmetry is spoiled by the renormalization. In the continuum 
the difference above is infinite, but we may use this formula on the finite lattice, 
giving 

k — 1 ■ ^ 

(23) u;iA;X) = Sxdk{A) = SxT^tv ^'^^' TK 

— 1 

Inserting T = ^4-+ Sx4- — [-4-^ + 1^^ (23) we obtain once more the 

form (19) of the anomaly. 

In the case of periodic boundary conditions in the continuum version (or more 
generally, with an enough rapid decrease of the gauge fields at |a;| oo) we can 
prove that (23) has a continuum limit as the momentum space lattice size is in- 
creased (the cut-off \p\ < A is removed). This follows from a simple Holder inequal- 
ity argument: If W is any operator such that \W\"' is a trace-class operator, that is 
W e Ln, and Wi is a sequence of operators in L„ converging to W with respect to 
the Ln norm then limtr(Wt)'^ = tiW^. Now W = T and the W^'s are z x z lattice 
approximations to the continuum operator T, lim \ W — Wj|„ = 0. 

The operators under the trace map in (23) become conditionally trace-class for 
k > 2n, [LMl], and the trace of the infinite dimensional matrices (continuum limit) 
is given by a local formula, the eq. (21) in the case 2n = 2. 

3. HAMILTONIAN FORMULATION: ANOMALY OF THE CURRENT AL- 

/^TTiTDD A 
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The same finite lattice approximation can be used also in the real time hamil- 
tonian formulation for fermions in background gauge fields. The hamilton operator 
Da in 2n — 1 dimensions is defined exactly in the same way as the 2n Dirac operator 
earlier, 

(24) Da = lQlk{Pk + Ak) = ak{pk + Ak), 

sum over k = 1, 2, . . . , 2n — 1. The dimension of the space lattice is now K = 
2^ ■ N ■ (2g + 1)^"^"^, in case of Dirac fermions; for Weyl fermions the hamiltonian 
above must be multiplied by the chiral projection P+ and the dimension of the 
projected subspace is K = 2"-^ -N ■{2q+ 1)^"-^ 

We shall consider the gauge currents for massless fermions. The left and right 
components of fermions decouple and we may restrict to (left-handed) Weyl fermions. 
The current algebra for Dirac fermions becomes just the direct sum of left and right 
current algebras. 

In continuum the current algebra is anomalous. There are Schwinger terms which 
in the case 2n — 1 = 1 can be written as 

(25) [p{X),p{Y)] = p{[X,Y]) + ^ [ trX{x)Y'{x)dx, 

27r Jm 

where the integral is over the one dimensional space M. When M is a unit circle 

(25) gives an affine Kac-Moody algebra. Here p{X) denotes the charge density 
integrated with a smooth Lie algebra valued test function X, 

p{X) = / pk{x)Xk{x)dx, 
Jm 

where A; is a Lie algebra index, k = 1,2, . . .,dimg. The trace under the integral 
sign in (25) refers to the representation of the gauge group acting on fermion com- 
ponents. In the case 2n — 1 = 3 one has, [Ml], [F-Sh], 

(26) [p(x), p(y)] = p([x, y]) + ^ ^ tTA[dx, dv], 

where the 3-form under the integral is defined as an exterior product of the 1-forms 
A, dX, dY. 

There are alternative, but equivalent, formulas for the Schwinger terms. Equiv- 
alent means again that the difference between the Schwinger terms is a gauge vari- 

— i- r ol A. \^\ r al A. v\ al a.W \/-\\ c ™„ .j-; — a A 
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and X, linear in the latter argument. In the case 2n — 1 = 1 one has in fact an 
exact formula, [L], 

(27) c{X,Y) = ^tre[e,X][e,Y], 

where e is the sign of the free 1-particle hamilton operator. In three space dimen- 
sions we have, [MR],[LM1], 

(28) c{A; X, Y) = itrc(e - F)[[e, X], [e, Y]], 

where now F is the sign of the 1-particle hamiltonian Da in the external gauge 
field A. The conditional trace is defined as trcK = |tr(ei^ + Ke). 

These formulas, and the corresponding formulas in higher dimensions, can be 
directly translated to the lattice. When defining the sign operators F, e one should 
be careful in order to have the right continuum limit. We have to split (somewhat 
artificially) the energy spectrum even in the finite case to positive and negative 
parts. As before we define the momentum components as Pi = -^^k with k = 
0, ±1, ±2, . . . , ±q for some positive integer q. The energy eigenvalues for the free 
hamiltonian become then E — ±[piPi)^/'^ and so the spectrum is symmetric around 
zero. 

In the case of periodic boundary conditions (in x space) in the continuum theory 
it is simple to prove that one obtains the right continuum limit as the size of the 
momentum space lattice is increased, that is, when the cut-ofi' \E\ < A is removed. 
This follows again from a standard Holder inequality argument. The trace in (28) is 
convergent, it is known that the operators [e, X], [e, Y] belong to the Schatten ideal 
L4 of operators T such that |T|^ is trace-class; furthermore, the diagonal blocks 
e(e — F) + (e — F)e oie — F (the only part of the operator which contributes to (28)) 
are Hilbert-Schmidt. All these operators can be approximated by finite-dimensional 
matrices (in the appropriate norms) and therefore the trace of the product is a 
limit of traces of finite-dimensional (cut-off) matrices. 

We have discussed above the abstract commutation relations of the current al- 
gebra but we have said nothing about an operator realization of the currents. In 
1-1-1 dimensions the theory is well understood; a physically acceptable realization 
is obtained using highest weight representations of affine algberas. An important 
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space. This construction has been generahzed to the 3 + 1 dimensional case (and 
the method works in higher dimensions) , [M2] . That representation is also suitable 
for a lattice approximation, as will be breafly explained below. 

The basic idea is to define a continuous family Ta of unitary conjugations in 
the one-particle fermionic Hilbert space such that the off-diagonal blocks (with 
respect to the energy polarization e) of currents are reduced such that the resulting 
(unitarily equivalent) Gauss law generators can be quantized by canonical methods. 
More precisely, we have to require that 

(29) [e, e{X; A)] is Hilbert- Schmidt , 
where X is any infinitesimal gauge transformation and 

(30) e{X;A) = Ta-\5x+X)Ta - Sx- 

Note that the modified Gauss law generators Gx = Ta~^GxTa, with Gx = Sx+X, 
automatically satisfy the same commutation relations as the generators X. The con- 
struction of the operators Ta is best understood through an (asymptotic) expansion 
in powers of the inverse momenta 1/^. In three space dimensions the first terms of 
the symbol are, [M2], 

(31) T^ = l_li[^,4]i + ... 
and the resulting gauge currents 

(32) e{X; A) = X+'-^[^, a^]d^X^ + ..., 

where we have used the hermitean Pauli matrices cr^ as the 3-space Dirac matrices, 
7^ = o'lj.Piu.- In all formulas it is implicitely assumed that an infrared regularization 
l/\p\ ^ l/(|p| + A) is performed. The terms which are of order strictly lower than 
— 1 in momenta are Hilbert- Schmidt and therefore not critical for the current renor- 
malization. With the renormalized current operators in hand, one can compute the 
quantum commutation relations in a straight-forward way. The resulting Schwinger 
term was calculated in [M2] and found to be equivalent with (26). 

These formulas can be translated immediately to the lattice. One replaces the 
derivatives —d^X by [p^, X] and the momentum symbols in the PSDO's become 
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